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Large deviations for directed 
percolation on a thin rectangle. 
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Following the recent investigations of Baik and Suidan in [I] and Bodineau and Martin in [5], we prove 
p> | large deviation properties for a last-passage percolation model in 1? + whose paths are close to the axis. The 

results are mainly obtained when the random weights are Gaussian or have a finite moment-generating 
function and rely, as in pQ and [5], on an embedding in Brownian paths and the KMT approximation. 
The study of the subexponential case completes the exposition. 



^ \ Introduction 

Random matrix theory has developed extensively in the last several decades following the pio- 
neering results by E. Wigner in the fifties. Gaussian models attracted a lot of attention, among 
them the Gaussian Unitary Ensemble (GUE). In this example, the knowledge of the joint dis- 
tribution of the eigenvalues allowed for a rather complete understanding of both their global 
and local behaviors. In particular, the limiting behavior of the largest eigenvalue gave rise to 
the famous Tracy- Widom distribution |15 | \2\\ l3T]. Several random growth models, such as the 
longest increasing subsequence and the corner growth models, have been shown to develop a 
similar behavior relying on a common determinantal structure |16| 118] , In particular, the last- 
passage percolation or so-called corner growth model (see below for a precise description) has 
been deeply studied by Johansson in |16| . For geometric or exponential random variables (the 
only cases leading to a determinantal description), Johansson established both fluctuations and 
large deviation asymptotics similar to the ones for the GUE random matrix model. Following 
the recent investigations by Baik and Suidan flj and Bodineau and Martin [5] at the level of 
fluctuations, the present paper deals with large deviations for the random growth model for more 
general random variables but on rectangles such that one side is asymptotically negligible with 
respect to the other at a given rate. The main results concern Gaussian random variables and 
random weights having a finite moment-generating function. Somewhat surprisingly, the rate 
may be shown to be larger than the one for the fluctuations. The comparison method used in 
this work is basically inspired from [1] and [5] and relies similarly on an embedding in Brownian 
paths. 

Recall first the basic corner growth model under study. It can be described as directed paths 
in the lattice going from (1,1) to (N,k) £ where only up and right steps are allowed. 
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More precisely, denoting by Tl(N, k) the set of all such paths, a path tt G II (iV, k) is called an 
up/right path and is defined as a collection of sites {{ii, ji)}[=[ k ~ 1 satisfying (i\,ii) = (1,1), 
(i N+k -i,jN+k-i) = (N,k) and - (U,ji) is either (1,0) or (0,1). The main random 

variable under consideration is the last-passage time defined by 



G(N,k) = max j V xj j) \ 
wen{N,k) L.^ J 



(jj')eTr 

where the X- '^'s are i.i.d. random variables. As an alternate description, set U(N,k) as the 
subset of B^ +1 given by 

U(N,k) = {u = {u ,ui, . . . ,u k ) G M+ +1 ;0 = u < ui < . . . < u k = N}. 

Then 

G(N,k) = sup {^[S^-S^ 



ueU(N,k) 

where Sm = with the convention S_"i = = 0. This follows from the fact that 

U(N, k) = U(N, k) when u G Z+ +1 . Actually, every u G t/(A, fe) n maps to a unique path 

7r G II(A r , k) whose i th up-jump occurs on On the other hand, each path 7r is characterized 
by its up-step sites. It will be more appropriate to adopt the second form of G(N, k) in order to 
compare it later with the Brownian last-passage percolation model. 



1 Results 

Having introduced the model, recall briefly some of its properties established by Johansson in 
|16] in the particular case of geometric or exponential distributions. The key of the results in 
this case relies on the explicit description of the last-passage time distribution G(N,k). When 
weights are i.i.d. geometric with parameter q G (0, 1), we have 

nG{N , k) < t]= j- E n <fc-*,)»n(*« + £-*)/., 

heN k l<i<j<k t=l V * ' 

max{hi}<t+k— 1 

where N > k and Zj^k is the normalizing constant. Using results from logarithmic potential 
theory, Johansson described in [16J the large deviation behaviors of G(N, k) when k ~ N. Namely, 
he obtained that, for 7 > 1, there exist two functions i(e) and /(e) such that for any e > 0, 

lim llogP[G([ 7 A], N) > N(u>fr,q) + e)] = -i(e) 

TV— >oo iV 

and 

lim -LlogF[G([ryN],N) < N(oj( 7 , q) - e)] = -1(e). 
The functions l(x) and i(x) are positive for every x > 0. Furthermore, 

lim 1 E[G([ 7 N}, N)] = a,( 7 , q) := (1 + - 1. 

7V->oo Jy 1 — 



1 Results 



3 



Using the asymptotics of the Meixner orthogonal polynomial ensemble, Johansson |16] further 
established the fluctuations of G(N, k) at the Tracy- Widom GUE rate. He proved that for 7 > 1 
and s£R, 

-G([ 7 iV],A0-AM7,9) 



lim I 

N-+oc 



cr(7,g)iVi/3 



< s 



Ftw(s), 



where 



°{l,q) = 1/6 (^+Vg) 2/3 (l + V79) 2/3 



and -Ftw( s ) is the distribution function of the Tracy- Widom law (see |31|). Replacing geometric 
weights with exponential ones gives similar results since an exponential distribution can be seen 
as the limit of a rescaled geometric one. See [16J for the precise formulas. 

Recently, Bodineau and Martin [3J and Baik and Suidan [JJ studied the same model when paths 
are close to the axis, i.e. k = o(N a ) for some a < 1 but allowing more general distributions. 
The authors used a coupling with the Brownian trajectories through the following Brownian 
last-passage percolation. Letting (B^ ) r >i be a sequence of independent Brownian motions, set 



L(N,k)= sup \j2l B $- B{ uL))- 



It has been proved in [2j [TJl [23] that L(l, k) has the same distribution as the largest eigenvalue 
of a k x k rescaled GUE random matrix. As a consequence of the fluctuation result for the GUE 
model, it follows that 

k 1/6 [L(l,k) -2Vk] AFtw- 

Using this result and a comparison between the continuous model with Brownian paths and the 
discrete one with random weights, the authors of [lj and [5] deduced fluctuation properties of 
the corner growth model for rather general random variables. However, the embedding in the 
Brownian paths requires to restrict the paths on small rectangles. For example, in [3], the discrete 
and the continuous models were coupled using the Komlos-Major-Tusnady (KMT) approximation 
which couples random walks with Brownian motion. The authors proved that if the weights satisfy 
E|X^|P < 00 for some p > 2, setting fi = EX^ and a 2 = var(X^ ) ), then for all a < f (5 - i), 

G(N, [N a \)-N^-2aN 1 ^ d 

j— ^ > ^TW- 

criVa 6 

If the random variables have all moments, i.e. p = 00, then a is lower than 3/7. This is true 
when the weights are Gaussian or are bounded for example. It is not known how optimal this rate 
could be : the authors in [5J think that such a result might hold, for some independence reasons, 
when a < 3/4. However, they do not give a complete proof. In [1], the authors compared the 
discrete and continuous model via the Skorokhod embedding theorem in order to obtain almost 
the same results. Lately, Suidan in |30| produced another proof of the last theorem when the 
variables have a third moment. He compared two discrete directed percolation models using a 
theorem of Chatterjee [8]. The fluctuation properties for the first (with geometric distribution) 
lead him to similar ones for the second. 
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In this paper, we follow the comparison methods of [5] and [T] to establish large deviations limit 
theorems for directed percolation models on thin rectangles mainly for Gaussian distribution. 
An extension for random weights with finite moment-generating function aroud zero will be 
also proved. We rely similarly on the corresponding results for the Brownian percolation model. 
Namely, as a consequence of the GUE random matrix interpretation |14l 124] . for all e > 0, 



lim \ logP[L(l, k) > 2Vk(l + e)] = -J GUE (e) (1.1) 



>oo k 
and 

lim ^ logP[L(l, k) < 2Vk(l - e)] = -I G ue(s). (1.2) 
fc-s-oo k z 1 J 

The two functions Jgue(x) and Igue(x) are both positive for every positive x. Jgue can be 
computed explicitly (see [3]) as 



Jgue{z) = 4 / y/x(x + 2) dx. 
J o 



To the best of our knowledge, there is no explicit form for Ique- This function appears in the 
logarithmic potential theory and it represents physically the minimal potential energy of charges 
on a one-dimension conductor exposed to an external field (sec |26j). In this work, we do not 
need the explicit form of Ique- However, its continuity, proved at the end of Section [2l will be 
necessary for the proof. 

The following three theorems are the main results of this paper. Despite some similarity in their 
proofs, the second one requires more work. One surprising feature is that the rate a in the 
Gaussian case is less then one (compared to a < 3/7 for the fluctuation result). 
Throughout the article, k and N are two integers which depend on each other, i.e. N = N(k) = 
iVfc. We assume that k = o(Nk) and we let k goes to infinity. For simplicity, we write ./V instead 
of iVjfc throughout the proofs. 

Theorem 1.1. Assume that the variables (X^)fj =1 are i.i.d. standard normal random variables. 
Assume further that k = o( lo ^- )• Then, for all e > 0, 

lim \ logP[G(^ fcj fe) > 2^^(1 + 6)] = -Jgue(s). 

fc^oo k 

i 

On the left of the mean, we have for k = o(N£ ), 

lim i log¥[G(N k ,k) < 2 V / AW - e)] = -I G ue(s). 

In the second statement, we replace Gaussian variables with weights having finite exponential 
moments. Loosing the Gaussian assumption will complicate the coupling and reduce the size of 
the rectangles. We denote by X a random variable having the common law of the i.i.d. variables 
in the sequence (X^^ CK 



Theorem 1.2. Assume that the variables {Xf )f- =1 are i.i.d. random variables such that EX 
and EX 2 = 1. Assume further that there exit fiQ > such that for all \i < hq, 



Eexp(^|X|) < +oo. 



(1.3) 
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lim 1 logP[G(^ fcj fc) > 2^/^(1 = -J G[/B (e). 

1 

Similarly, if k = o(N£), for all e > 0, 

lim i logP[G(iV fc ,&) < 20V^(1 - e)] = -I GC/£ (e). 

/C — ^OO fx, 

The proof of Theorem 11.21 relies on the Komlos-Major-Tusnady approximation for the sums of 
i.i.d centered random variables with finite exponential moments, see |17j . The following theorem 
deals with a particular class of subexponential weights. We make use of the Skorokhod embedding 
theorem to obtain the following results. 

Theorem 1.3. Assume that the variables (X^ ; )^- =1 are i.i.d. random variables satisfying EX = 
and KX 2 = 1. Furthermore, assume that there exit fj, > and < 7 < 1 such that 

Eexp(//|X| 7 ) < +00. (1.4) 

If k = o(N£) with a < 2 ^+2 > then, for all e > 0, 

lim \ logP[G(N k ,k) > 2^/N^k{l + e)] = -J GUE (e). 
Similarly, if k = o(N£) with a < 5 ^_ 4 , for all e > 0, 

lim ^ logP[G(JV fcj fc) < 2^fN~kk{l-e)] = -I GU e{b). 

The results in Theorem 11.31 cover in particular the examples of Weibull and Levy distributions. 
Notice that in Theorem 11.31 we can take 7 > 1. However, the result is worthless because of 
Theorem 11.21 Actually, the KMT approximation is more efficient than the Skorokhod embedding 
theorem for i.i.d. random variables with finite exponential moments. 



The asymptotic results above hold also for geometric and exponential weights. Moreover, the 
rectangle width in this case, has to be only negligible with respect to its length since the rate 
functions of the Laguerre ensemble converge to the GUE ones when k = o(N/ t ). The reader can 
see |16] and |19) for rigorous results. 

Non-asymptotic bounds for the preceding models can be deduced from the previous theorems 
proofs. The rectangle width for small deviations matches in this case the fluctuation results. For 
this we use analogous deviation inequalities to the right of the mean obtained for the largest 
eigenvalue of the GUE , see |19| . To the left of the mean, we use recent deviation results for the 
largest eigenvalue of the GUE obtained by Ledoux and Rider |20j. 

Theorem 1.4. Assume that the variables (X^)fj =1 are i.i.d. standard normal random variables, 
and that k = N£ with a < I. Then, there exists a positive constant C a depending only on a such 
that, for all < e < 1, 

, / fee 3 / 2 \ 

P[G(N, k) > 2y/Wk{\ + e)] < C a exp ( - — — ) . (1.5) 
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On the left of the mean, we have 

i 2 3 

P[G(N,k) <2VNk(l-e)] <C a exp(-— (1.6) 

For random weights with finite moment-generating function, we have similar results. 

Theorem 1.5. Assume that the variables (X^)f 3 j =1 are i.i.d. random variables such that EX = 
and EX 2 = 1. Assume further that there exit fj,Q > such that for all pL < /jlq, 

Eexp(^|X|) < +oo. (1.7) 

. If k = N£ with a < S, i/ien i/iere exists a positive constant C a depending a and the distribution 
of X such that, for all < e < 1, 

ke 3 ' 2 \ 



F[G(N k ,k)>2y^k(l + e)] <C a exp(^-— J. (1.8) 
Similarly, if k = with a < \, then 



k 2 e 3 



F[G(N k ,k) < - e)] < C a exp ^ - —J . (1.9) 

Large deviation inequalities for large e > 1 hold for the optimal rate a = 1 using some concen- 
tration arguments. We refer to |19| for the results and the proofs. 

As the reader can notice, if X satisfies the condition (|1.4I) . the same exponential inequalities as 
in Theorem 1 1 . 51 can be obtained. In this case, the positive constant will depend on a, [i and 7. 
A smaller a will be also necessary. The precise calculations are left to the reader. 
Let us notice here that if the weights have a finite p-th moment with p > 2, we still obtain 
the same deviation results but for k = o(log(N p ~ 1 )). In this case, we use again the KMT 
approximation for random weights with finite p-th. moment, see [171 127] . 

We strongly believe that the preceding results hold on a wider rectangle. However, the method 
used here does not allow us to improve the rectangle width. Theorem II. II will be proved in Section 
[5] while Theorem 11.21 will be proved in Section |3J The case of the particular subexponential 
weights given in Theorem 11.31 will be discussed in section 21 Theorem 11.41 and Theorem 11.51 will 
be addressed in Section [5] on the basis of the preceding results and proofs. 



2 Proof of Theorem 11.11 

We assume throughout the rest of the paper and without loss of generality that k and iV are two 
positive integers and we write instead of Nf. and N(k). As claimed before, to prove Theorem 
II. 1| we compare G(N, k) and L(N, k). To do so, let for any e > 0, 

A = {G(N, k) > 2VNk{l + e)} 

and 

B = {\G(N,k)-L(N,k)\ > 2^Nk{e-ei)} 
where < e\ < e. Clearly, 

F[A] < F[L(N, k) > 2v / Afc(l + £1)] + F[B] (2.1) 



2 Proof of Theorem [7771 
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and 

F[A] >F[L(N,k) >2\/iVfc(l + 2e-ei)] -F[B]. 
Moreover, for every rj > 0, 

F[L(N, k) > 2VNk(l + rj)] = F[L(1, k) > 2y/k(l + rj)] 

d 



(2.2) 



(2.3) 



as a consequence of the Brownian scaling y NL(l,k) = L(N,k). To evaluate F[B], we couple 
G(N, k) and L(N, k) by letting X^ j) = B^-B^ for all i,j>l so that the sequence (X i (j) )~ =1 
is i.i.d. with standard normal distribution. When comparing G(N,k) and L(N,k), it is obvious 
that most of the variables will vanish. More precisely, repeating the computation done by Bodi- 
neau and Martin in Section 2 of [5], we get, by letting B_± = 0, 

\G(N,k)-L(N,k)\ 



sup , I^M w K-iJ-i 



sup 

u'eu(N,k) r= i 1 Ur 



r) _ B (r) 



< sup Yl \\ 



n(r) _ R (r) 
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(r) 



13 



(r) 



+ \b\ t \-b^ 

\U r \ u r 



+ 



|u r _lj-l L"r-lJ-l 



< 



2^( max |5«-^|)+2^( sup |bW-^|) 



Let us try to find an exponentially decreasing upper bound for each term of ([2 

k 



(2.4) 

. For this, set 



and 



Then 



r=l 



( max 


S r>-Br> 


V 8=1,. ..,JV 





Z k = 2^( sup B^-bP). 

1 — * \ 0<s,t<N I 

r=l v |7-t|<a 

P[S] < P[E fc + ^ > 2\/jVfe(e-ei)] 

< F[Y h >VNk(e-e 1 )]+¥[Z k >VWk(e-e 1 )]. 



(2.5) 
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In view of the Gaussian hypothesis and the coupling, Y k = 0. Applying the Markov inequality 
gives for any A > 0, 

F[Z k > VWk(e-ei)] < E[exp (AZf)] . exp ( - A(e - £i) 2 AA;) 

exp(2Afc( sup [BP - B { t 1) \Y)] k .exp(- X(e-ei) 2 Nk) 

V V 0<s,t<N ' /J 

|s— 1|<2 

J °°4Afct.exp(2AA:t 2 )J 



However, 



< E 



< 



sup \Bf ) -B [ i ) \>t 

0<s,t<N 
\s—t\<2 



dt 



exp(X(e-ei) 2 Nk) 



N-3 



sup \BP -BP\ >t] < X) P [ SU P B t~ .J5f,„ s * ^* 



0<s,t<N 
|s-t|<2 



i=0 



-i<t<i+3 i<t<i+3 



(2.6) 



< NF 



sup \B t \ > t/2 

0<t<3 



By the Brownian motion reflection principle (see for example |25|). snp 0<t<a Bt = \B a \. Thus, 

sup \BW - b[ 1] \ >t] < 4Nf\b 3 > t/2 



0<s,t<N 
\s— 1|<2 



< CiAexpf - — 



(2.7) 



where C\ is a numerical positive constant. Now, we insert (|2.7p in the integral in (|2.6p and we 
choose A = ^ were c is a positive constant smaller than Then we get 



\Z k > VNk(e - £i)] < C 2 exp ( ° ) . exp 



C 2 



(e-ei) 2 N 
C 2 



(2.{ 



The last bound leads to the condition k = o(N/\ogN). Combining ([2TT]) . (l2~2l) . ([231) and ([2~8]) 
then leads to 

P[A] < P [L(l, fc) > 2Vk(l + ei)]+ C 2 exp ( - (g ~ £l) ^ - kl ° gN ) (2.9) 



and 



PL4] >P[L(l,fc) > 2^(1 + 26 -£i)] -C 2 exp( 



(e-e^N -klogN 

Co 



(2.10) 



2 Proof of Theorem ll.il 
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Dividing fl2j| by e ^ kJauE ^ and (f2~T0]) by e^-W^e-ei^ ta k ing t h eir logarithm and then 
dividing the results by k, we get for a < ^, 

1, / F\A] \ 1 fF\L(l,k) > 2Vk(l + ei )1 A . 

-log = / J . , < - log 1 V ' 7 - , \ ;J +gfc £i,£ (2.11 
k \ e - kJ GUE(ei) J ~ k \ e -kJ G u E (ei) > J v ; 

and 

1 / F[A] \ 1 fP[L(l,k) > 2y/k(l + 2e-£i) 



1, ( \ ^ 1 1 ^ P[L(l,fc) > 2y/fc(l + 2£-£i)] , A 



where 



v „ ( (e-ei) 2 N , AogN 

ffjfe(ei,e) = C 2 exp ^ ^ Jgue(si) 

and 

^(ei.e) = C 2 exp( - ^ + + J GC / £ (2e - £l ) 

are two positive functions. Moreover, for k large enough, gk(ei,e) and g((ei,e) are negligible 
with respect to e~ vk for every r/ > since & = o(N/logN). Thus, using (|l.ip . a straightforward 
computation shows that the right-hand sides of (|2.1ip and ()2. 12[) both converge to zero when 
k — > 00. In other words, for k = o(N/ log N) and E\ < e, 

limsupilogP[A] < -J GUE (ei) (2.13) 

k— >oo ™ 

and 

liminf ^logP[A] > -J GUE {2e-e 1 ). (2.14) 

Finally, notice that Jgue(z) is a continuous function of e > 0. It therefore follows from (|2,13p 
and (I2.14P that for every e > 0, 

lim jlogF[G(N,k)>2VNk(l + 8)}=-J GUE (e). 

k— >oo K 

The proof of the leftmost charge formula is similar. Set now, for all e > and E\ < £, 

E = {G(N, k) < 2\/jVfc(l - e)}. 
By the same arguments as before, we get 

F[E] < F[L(N, k) < 2VWk(l - e x )] + F[B] (2.15) 

and 

F[E] > F[L(N,k) < 2VWk(l -2e + ei)] -F[B}. (2.16) 

Furthermore, by (|1.2p . 

1, /F\L(l,k)<2Vk(l-ei)]\ 
lim -r log 1 K — '-. \ ^ = 2.17 

fc^oo k 2 V e ~ k T GUE(ei) I 
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and 



lim 1 log ( ^(1,^ 2^(1-2^)] ) _ 

fc->oo k 2 V e -fc 2 I G <7i3(2e-ei) ' 



Using the same upper bound on ¥[B] and combining (|2,15p . (I2.16p . (|2.17p and (I2,18p . one can 
easily deduce that, for k = o(N 1 ^ 2 ) and E\ < e, 

lim sup log ¥[E] < -lGUE(ei) 

k— >oo ™ 

and 

liminf rjlogP^ > -I GUE {2e-e l ). 

A;— >oo fe z 

At this stage, let us assume that Igue{z) is a continuous function of e. Then, for k = o{N 1 / 2 ), 
lim 



lim -L logP[G(JV,fc) < 2x/]Vfc(l- £ )] = -I GUE (e), 



which is the result. 

We are left with the proof of the continuity of Igue{^)- Set oo,t]), the set of all proba- 

bility measures on (— oo,t] when t G R. For a given distribution /i G oo,i]), define the 

corresponding potential energy, as in |26j . by 



I »(t) = 2 l x 2 dfi(x)- / log \x - y\dfi{x)d/i(y). 



oo ■/ — oo 



The minimal energy 



I(t) = inf I„(i) 

fj,eM{(-oo,t]) 



precisely allows us to compute the rate function Igue via the formula Igue(^) = I(l—e)—I(oo). 
The last equality could be found in [llj. For t > 1, I(t) is a constant function, the extremal 
measure is the so-called semi-circular law supported on [—1, 1] and the energy I(t) = log (2) +3/4, 
(cf. [H [26]). For each t G R, there is a unique measure vt G ,M((— oo,i]), with no mass point, 
achieving the infimum (cf. [26]). Furthermore, vt is compactly supported and the corresponding 
energy is finite. Since I(t) is an infimum and a non-increasing function of t, for any r\ > 0, 



Iu t (t-V) 
u?((-oo,t-rj\) 



i(t) < i(t - V )< . 2 ,/*:_ , _n • ( 2 - 19 ) 



It is obvious that the right-hand side of (I2.19P converges to I(t) when r\ converges to zero. This 
proves the left-continuity of l(t). 

To show the right-continuity, notice that by a simple change of variable, 

/(*)= inf Mt) 
neM((-oo,t+ij\) 

where 

/t+r) rt+rj rt+rj 

{x- V ) 2 dfi{x)- / log\x - y\dfi(x)dfi(y). 

-OO J —oo J —OO 



3 Exponential-tailed distribution and the KMT approximation 
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Consequently, 

I(t)-I(t + rj) <II +V (t)-I(t + V ) 

rt+v 

<2r] 2 + Ar] \x\ du t+v (x). (2.20) 

J — oo 

For \x\ > \y\, we have log \x — y\ < log \2x\. Moreover, there is a positive constant C3 such that 
\x\ < Cs(2x 2 — 2 log \2x\). In view of (I2.20p . a straightforward calculation leads to 

I(t)-I(t + rj) <2r] 2 + 4 V C 3 I(t). (2.21) 

Since I(t) is finite, the right-hand side of (|2.21[) converges to zero when 77 — > 0. Thus, the 
continuity of I(t) is proved, and that of Igue{^) as well. The proof of Theorem 11.11 is now 
complete. 

□ 



3 Exponential-tailed distribution and the KMT approximation 

In this section, we replace the standard normal variables with weights having finite moment- 
generating function around zero. When comparing G(N, k) to the Brownian last-passage per- 
colation model, Yk will not vanish as in the Gaussian case where the coupling was "perfect". 
Actually, we couple easily a partial sum of i.i.d. standard normal random weights with a Bro- 
wian motion. However, loosing this assumption will complicate the task and a new coupling is 
then required. Following [5], we make use of the KMT approximation : a powerful tool to couple 
a partial sum of i.i.d. random variables and a Wiener process, both constructed on the same pro- 
bability space. The KMT approximation, also called the invariance principle, was first introduced 
in 1975 by Komlos, Major and Tusnady in the famous work [17] , The basic version deals with 
a partial sum of i.i.d. random variables reconstructed in a way to be "close" to another partial 
sum of i.i.d. standard normal random variables. Later versions of this strong approximation do 
not require a common distribution, see |27| . The readers can also see |9j for a complete survey. 
Let (Xi)i>i be a sequence of independent random variables and denote by Sn the corresponding 
partial sum. Let (Bt)t>o be a Brownian motion built on the same probability space. The following 
theorem is an immediate consequence of Theorem 1 in [T7] . 

Theorem 3.1 (Komlos-Major- Tusnady). Assume that KXi = and KX 2 = 1. Assume further 
that there exit fiQ > such that for all /i < 

Eexp (n\Xi\) < +00. 

Then for every N > 1, the sequence (Xj)j>i and the Brownian motion (Bt)t>o can be constructed 
in such a way that for all x > 0, 

¥[ . max \S l - B,\ > GlogiV + x] < Cexp(-6x). 

i=l : ...,N 

The positive constants Q, C and 9 depend only on the distribution of X\ and can be taken as 
large as desired by choosing large enough. 
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According to the notation and the steps of Section [2j recall that 

\G(N,k)-L(N,k)\<Y k + Z k . 



We already have Z k in fj2.8j) and we want Y k to be as small as possible. To this, we construct the 
sequence (X^)f°j =1 and the independent Brownian motions (B^ )t>o i n the sense of Theorem 
3.11 By the Markov inequality we have for all e > 0, E\ < £ and A > 0, 



} [Y k > v / M(e-ei)] < E[exp(AYfc)].exp(- A(e-ei)v / iVfc) 

exp ( 2 A max 

. V 1=1,...,, 

J °°2A.exp(2At). 



< E 



< 



.exp ( - A(e-ei)VFfc) 
fe 



maxi=i jv 



> t 



dt 



exp (A(e - ei)VM) 



In order to apply Theorem 13.11 we make the simple variable change t = s — G log N and we 
choose A < 9/2. Therefore, there exist two positive constant C4 and C4 such that, 



>[Y k > ^Nk(e - ei)] < C 4 exp 



(e-e^VNk- klogN 



Now, putting (pTTj) . (J2T2J), (gJJ and (|3T]) together gives 

P[A| <F\L(N,k) > 2v / iVA?(l +£i)l +C 5 exp 
(e-ei)VNk- klogN 



( e -£i) 2 iV-HogAf 



(3.1) 



(3.2) 



C5 exp 



Co 



and 



F[A] >P L(N, k) > 2VNk(l + 2e- £1) 



C5 exp 



(£-£l) 2 iV- klogN 



C5 exp 



(e-e-fiVNk- klogN 

C~5 



(3.3) 



On the left of the mean, we have 

P[E] <p[l,(iV,fc) < 2VNk(l-e 1 )] +C 5 exp 
(e-ei)VWk- klogN 



(e - ei) 2 N -klogN 



(3.4) 



+ C 5 exp 



C 5 



and 



¥>[E] >P L{N, k) < 2VNk(l -2e + 



C5 exp 



(£-£i) 2 iV- klogN 



— C5 exp 



(e-ei)\/jVfc-fclogJV 
C 5 



(3.5) 
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Proceeding like in Section we divide (JUS) by e - kjGUE ^ , §33$) by e -* J GUE(&-ei) ^ (jg^J) by 
e -fe Igue(£i) anc [ (|3,5p by e _fc fct/M2e-£i)_ ^ bandle the remaining parts when k — > oo, we take 
k negligible with respect to N. On the right of the mean, we need k = °( /i ^j\As ) an d on the 
left, we take k = o(iV 1 / 3 ). Finally we conclude as in section [2] using the continuity of J que and 
Igue- 



4 Subexponential weights and the Skorokhod embedding 

In this section, we consider a particular category of subexponential weights verifying E exp (/i|X| 7 ) 
< oo for some fi > and 7 6 (0, 1). We say that these variables have a Weibull-like tails because 
of the similarity with the right tail of the Weibull distribution with a shape parameter lower than 
1. Such weights are considered to be heavy-tailed and then, do not satisfy the Cramer condition. 
However, G(N, k) still satisfies an LDP principle on a very thin rectangle. For subexponential 
weights with finite p-th moment for p > 2, we need k to be at least smaller than log (iV p_1 ) in 
order to prove deviation results. Computations in this case are left to the reader. The precise 
definition of a subexponential distribution and large deviations for a partial sum of i.i.d. such 
weights can be found in [22J. 

We do not know a KMT strong approximation version for random variables satisfying the mo- 
ment condition above. We follow [1] by using the Skorokhod embedding theorem instead. The 
Skorokhod embedding theorem is another tool to couple a sum of i.i.d. random variables with a 
Brownian motion, see (cf. [6J [29J [23]) f° r more details concerning this theorem. 

Theorem 4.1 (Skorokhod). Let (Bf)t>o be a standard one- dimensional Brownian motion and 
X a real valued random variable satisfying EX = and EX 2 = 1. Then, there is a stopping time 
T for the Brownian motion such that B? = X and ET = 1. 

An immediate consequence of this theorem allows to embed sums of real independent random 
variables into the Brownian motion. Applying the strong Markov property to the Brownian 
motion, Theorem 14.11 yields the following classical corollary. 

Corollary 4.2. Let X\, X2, ■ ■ ■ , Xjy, . . . be i.i.d. satisfying EX\ = 0, EX 2 = 1 and set Sn = 
X\ + X2 + • • • + Xn> N > 1. There is a sequence of i.i.d. stopping times tq = 0,T\, . . . , rjy, . . . 
such that 

Sn = B Tl ^ — |_ TJV 

and (.B T1 _| \-t n+1 — B T1 ^ \- tn )n>o is a sequence of i.i.d. random variables having the same 

distribution as X\. 

In our context, an application of the last corollary allows to claim that there exists i.i.d. stopping 

(r) d (r) 

times for the Brownian motion tq = 0, t±, . . . ,tn, . . . such that Eti = 1 and = j_ T . for 

i > 1 and r > 1. Consequently, choose 

y( r ) — dO) _ rM 

A % — n T 1 + -+T i i3 r 1 + --- + Ti_l! 



in order to have 



S i = B k+-+ri a - s - 
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Thus, for any t > 0, 



max \sP -sPl > t 



max B 



(i) 



Ml 



Li=l AT 



i=i 



> v 



+ 



max \Bill +T . - B\ l) \ > t; max I V (r, - 1) 



!=1 



Hence 



max Isf 5 -fif^l >t 
i=l,...,N 1 1 



< 



+P 



sup 

0<s,t<N 



max 
t=l,...JV 



BCD _ B W 



> t 



(4.1) 



We evaluate each term of (|4.ip separately. First, 



sup 

0<s,t<N 
\s—t\<tP 



> t 



N-tf> 



< 



V P [ sup S t - inf B t > t 

L »<i<i+f9+l i<t<i+t?+l 



i=Q 



< iVP 



sup |5 f | > i/2 



Applying the reflection principle as in Section [2j we get 



sup 

- a<s,t<N 

s-t\<tP 



BP - B® 



> t 



< 4NF 



BtP+i > t/2 



< AN exp 



(4.2) 



To find an upper bound for the second term on the right-hand side of (|4.ip . we need a connection 
between the weight moments and those of the stopping times obtained by the Skorokhod embed- 
ding. Furthermore, we need to control the sum of the independent stopping times to reach an 
exponentially decaying inequality. When the weights are bounded for exemple, we can construct 
a stopping time with finite exponential moments. The sum is then controlled by the Bernstein 
inequality. 
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However, when X only satisfies f)l -4[> . the Skorokhod stopping time does not necessarily have a 
finite exponential moment and thus the Bernstein inequality can not be applied. For example, in 
[10] . Davis found the best universal constant connecting the stopping time moments to those of 
the stopped Brownian motion. More precisely, if (Bt)t>o 1S a Brownian motion and r is a stopping 
time, then there is a universal constant a p such that, when 1 < p < oo and Et p / 2 < +oo, 

a p Er p/2 < K\B T \ P . (4.3) 

Moreover, Davis proved that the best constant for p = 2n (n £ N*) is z\^ n which is the 
smallest positive zero of the Hermite polynomial of order In. In [7], this constant is shown to be 
0((2n)~ n )). So unless X = B T is a bounded variable, r can not have finite exponential moments. 
The constant above is universal but it could be sharpened for some particular stopping times. For 
example, considering the stopping time of the Skorokhod representation [6, 29j, Sawyer improved 
the constant a p and established, in [28], the following inequality. 

Theorem 4.3. Let X be a centered random variable such that 

Eexp(^|X| 7 ) < +oo 

for some 7 > and u > 0, and let r be the corresponding stopping time of the Skorokhod 
representation. Set 6 = an d v = u 18 ■ Then, 

Eexp(z^T e ) < $ 7 Eexp (^|X| 7 ), 
for some positive constant <3? 7 depending only on 7. 

Note that a similar exponential bound may be obtained from (|4.3p . However the cost is a worse 
constant \i. 

Under the assumption (|1.4p and in view of Theorem 14. 3| the Bernstein inequality can not be 
applied to the sum of the independent stopping times because = 7/(2 + 7) < 1. To avoid this 
obstacle, we introduce the Fuk-Nagaev inequality [13] which requires less restrictive assumptions. 

Theorem 4.4 (Fuk-Nagaev). Let Xi,--- ,X/v be a sequence of real i.i.d. random variables sa- 
tisfying KXi = and KXf = a 2 . Then, for all x > and y > 0, 

„2 



max 

Ki<N 



1=1 



£X,|>X < W P|.Y,|>, + 2exp(- — £— (44) 



We refer to |12| for more detail on this inequality. Recall now the second term of the right hand 
side of (|4.ip and consider the stopping times of the Skorokhod representation. Choosing x = t^ 
and y = t s in Theorem 14.41 and applying Markov inequality to the first term of the right hand 
side of (|4.4|) . one has 

I ' , , v t es \ / t 2/3 \ 

35, £<* - >>M s C « N ■* ( - or) + c ° exp ( - owot) 1151 



1=1 



Now, we apply Markov inequality to as in Section [3] and we get for some A > andO < 77 < 1, 
P[Y k > VM(e - ei)] < E[exp(Ay fc ' 7 )].exp(-A(e-ei) r '(iVA;) ,? / 2 ) 



< 



/ 1 °°4Aexp(2At , ')l 



maxj = i v .. iA r 



<?« _ bP 



> t 



dt 



exp (A(e -e 1 )n(NkYl 2 ) 



(4.6) 
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Inserting (|4.ip and (j4.5|) in (j4.6|) and choosing A very small yield the following constraints on rj, 
P, S and 9. 

7] < 2-/3 
j] < @ — 6 
rj < 95. 

Straightforward computations lead us to choose /3 = 3^+2 an< ^ ^ = 37+2 smce ^ = ^+2- Conse- 
quently, we obtain r\ < 3 ^ 2 • ^° S e t large deviation asymptotic formulas on the right and the left 
of the mean, we respectively need a < and a < j^—- This completes the proof of Theorem 
IP 

□ 



5 Small and large deviations inequalities 

Non-asymptotic bound on the right and the left of the mean is an immediate consequence of 
the corresponding bound for the GUE and the arguments developed in Sections 2 and 3. In 
particular, we use (cf. [19]) that there exists a positive constant C7 such that, for any e > 0, 

P[L(1, k) > 2Vk{l + e)] < exp ( - U GUE (e)) < C 7 exp ( - fcmax gV 3/2 ) ^) . (5.!) 

On the left of the mean, deviation inequalities for the largest eigenvalue of the GUE for a given 
k are quite more complicated to prove. Ledoux and Rider obtained in a recent paper, [20J, that 
the leftmost charge of the largest eigenvalue of a large set of random matrices behaves like the 
left tail of the corresponding Tracy- Widom law. More precisely, they get for all < e < 1, 

P[L(1, k) < 2Vk(l - e)) < C 7 exp ( - ^\ . (5.2) 

As we mentioned before, when e > 1, we have Gaussian behavior for both left and right tails. This 
follows from concentration arguments dealing with Lipschitz functions of independent standard 
normal variables. Once more, two cases will be tackled : Standard normal weights and finite 
exponential moments ones. 



5.1 Standard normal variables 

Following the proof of Theorem 11.11 in Section [21 choose e\ = |. Then, combining (|2.9p and (|5.ip . 

for any e > 0, 

mr „ ( fcmax(e 3//2 ,£ 2 )\ ^ / e 2 N — klogN 
F[A] < C 8 exp ( ^ '—L + C 8 exp ' 



C% J V Cg 

where Cg > 0. In order to reach (|1.8p when < e < 1, we need a positive constant C(a) > Cg, 
depending only on a, such that 

/ ke 3 / 2 \ ( e 2 N- klogN ^ 
C a exp y - j > C 8 exp ^ J . (5.3) 



5 Small and large deviations inequalities 
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Taking the logarithm of f|5 . 3[> . C a has to satisfy 
However, since ¥[A] < 1, e has to satisfy 

fe 3/2 > 1. (5.5) 
Combining now f|5.4|) and (|5.5p . we finally get that C a has to satisfy 

C 8 1 N 1 '^ -AT" log AT 
!og 77- + r r °- ^ 5 - 6 ) 

Hence C Q exists and satisfies <j5.6|) only if a < |. In that case, we make the reverse computation 
to conclude that 

' [G(N, k) > 2v / iVfc(l + e)] < 2C a exp 



We make the same computations for the left-tail upper bound. Here, C a has to satisfy 

C 8 , 2 2 / iV lQgiV E\ 

which finally gives 

C 8 1 iV 1_2Q -iV Q/3 logiV 
!og 77- + 7^ 7; < °- 

This proves Theorem ll.4i 

□ 



5.2 Finite moment-generating function case 

Choosing E\ = | in (|3.2p and taking into consideration (|5.5p . the inequalities ([5.ip and (I5.2p 
imply that there exists a positive constant C9 depending on a and the distribution of X such 
that, for all e > 0, 

fe 3 / 2 \/ / JV2-f - iV Q log iV - iV c 



PM <C 9 exp ( - _) (l + exp ( ^ ) 

N l-4a/3 _ N a log N _ jya 



+ exp 



Co 



and 



F[E] <Cg exp - — - 1 + exp + exp 



Cg J \ V Cg J V Cg 

This means that we have a right-tail bound for a < 3/7 and a left-tail bound for a < 1/3. The 
proof is complete and thus Theorem 11.51 is proved. 

□ 
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